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Abstract 

We apply the ADM 3+1 formalism to derive the general relativistic magnetohydrodynamic equations 
for cold plasma in spatially flat Schwarzschild metric. Respective perturbed equations are linearized for 
non-magnetized and magnetized plasmas both in non-rotating and rotating backgrounds. These are then 
Fourier analyzed and the corresponding dispersion relations are obtained. These relations are discussed for 
the existence of waves with positive angular frequency in the region near the horizon. Our results support the 
[ fact that no information can be extracted from the Schwarzschild black hole. It is concluded that negative 

<sJ 'i phase velocity propagates in the rotating background whether the black hole is rotating or non-rotating. 
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1 Introduction 



General Relativity (GR) serves both as a theory of space and time as well as a theory of gravitation. This 
theory has tremendous philosophical implications and has given rise to exotic new physical concepts like black 

■ holes and dark matter. All the predictions and applications by different scientists make the theory a spec- 
I tacularly successful one. The theory of GR is not complete unless it includes the hard reality of experiments 

■ and of astronomical observations. To explain the phenomena observed by astronomers in physical terms we 
I have astrophysical relativity. This serves as a tool to gauge the properties of large scale structures for which 
I gravitation plays a significant role. 

It is believed that the compact objects have strong gravitational fields near their surfaces [1]. The study of 
. ■ general relativistic effects on electromagnetic processes which take place in the vicinity of compact objects is of 

r> I theoretical interest in its own right. However, it is also important to interpret and understand the still puzzling 

. signals we receive from these astronomical objects, including neutron stars and black holes. All the massive stars, 

neutron stars and black holes carry energy flux. Due to this flux, a relatively large magnetic field is produced. 
In order to understand the phenomena of the compact objects and many more like these, we have a developed 
theory of magnetized plasma called theory of magnetohydrodynamics (MHD). To study the magnetosphere of 
massive black holes, the strength of gravity demands the theory of general relativistic magnetohydrodynamics 
(GRMHD). GRMHD equations help us to study stationary configurations and dynamic evolution of conducting 
fluid in a magnetosphere. These equations (Maxwell's equations. Ohm's law, mass, momentum and energy 
conservation equations) are required to investigate various aspects of the interaction of relativistic gravity with 
plasma's magnetic field. The successful study of plasmas in the black hole environment is important because a 
successful study of waves will be of great value in aiding the observational identification of black hole candidates. 
It is well-known [2-4] that an isolated black hole cannot have an electromagnetic field unless it is endowed with 
a net electric charge. Since a collapsed object can have a very strong effect on an electromagnetic field, it is of 
interest to determine this effect when a black hole is placed in an external electromagnetic field using GRMHD 
equations. 

The 3+1 formalism (also called ADM formalism) was originally developed to study the quantization of grav- 
itational field by Arnowitt et al. [5] . Tarfulea [6] worked on numerical solutions of Einstein field equations by 
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applying constraint preserving boundary conditions. Thorne and Macdonald [7,8] extended the formulation to 
electromagnetic fields. Thorne ct al. [9] described its further applications for the black hole theory. The black 
hole theory was considered by Holcomb and Tajima [10], Holcomb [11] and Dettmann et al. [12] to investigate 
some properties of wave propagation in Friedmann universe. Khanna [13] derived the MHD equations describing 
the two component plasma theory of Kerr black hole in 3+1 split. Anton et al. [14] used 3+1 formalism to inves- 
tigate various test simulations and discussed magneto-rotational instability of accretion disks. Anile [15] worked 
on propagation and stability of relativistic shocks and relativistic simple waves in magneto-fluids. He also lin- 
earized the electromagnetic wave theory in cold relativistic plasma. Komissarov [16] discussed Blcndford-Znajec 
monopole solution using 3+1 formalism in black hole electrodynamics. In another paper [17], he formulated 
and solved time dependent, force free, degenerate electrodynamics as a hyperbolic system of conservations laws 
for pulsars. He also discussed the waves mode and one dimensional numerical scheme based on linear and exact 
Riemann solvers. Zhang [18] formulated the black hole theory for stationary symmetric GRMHD using 3 + 1 
formalism. The same author [19] applied this formulation to the linearized waves propagating in two dimensions 
using perfect GRMHD fluid for the cold plasma state in the vicinity of Kerr black hole magnetosphere. He 
investigated the behavior of perturbations in graphic form which described the response of the magnetosphere 
to oscillatory driving forces near the plasma-injection plane. Tsikarishvili et al. [20] considered 3+1 general 
relativistic hydrodynamical equations describing strongly magnetized collisionless plasma with an anisotropic 
pressure and obtained the equations of state in ultra-relativistic limit. In his recent paper, Mofiz [21] discussed 
gravitational waves in magnetized plasmas. 

The formalism for gravitational perturbations away from a Scwarzschild background was developed by Regge 
and Wheeler [22] . It was extended by Zerilli [23] who proved that the perturbations corresponding to a change 
in mass, the angular momentum and charge, of the Schwarzschild black hole are well-behaved. The decay of 
noil well-behaved perturbations was investigated by Price [24] . The quasi-static electric problem was solved by 
Hanni and Ruffini [25] who showed that the lines of force diverge at the horizon for the observer at infinity. 
Wald [26] derived the solution for electromagnetic field occurring when a stationary axisymmetric black hole 
is placed in a uniform magnetic field aligned along the symmetry axis of black hole. A linearized treatment 
of plasma waves for special relativistic formulation of Schwarzschild black holes was developed by Sakai and 
Kawata [27]. This treatment was extended to waves in general relativistic two component plasma in 3+1 ADM 
formalism, propagating in radial direction by Buzzi et al. [28] . They investigated the one dimensional radial 
propagation of transverse and longitudinal waves close to the Schwarzschild horizon. 

This paper is devoted to study the dynamical magnetosphere of the Schwarzschild black hole. The model 
spacetime is the Schwarzschild metric in Rindler coordinates which preserves the key features of the Schwarzschild 
geometry. This investigation supports the well-known fact that the information from the black hole cannot be 
extracted. We apply 3+1 form of the GRMHD equations to this spacetime by using the cold plasma state. 
Afterwards, we linearize the equations to study perturbations of the magnetosphere. These equations arc then 
Fourier analyzed to get the sets of ordinary differential equations which are comparatively easier to handle than 
the original GRMHD partial differential equations. The determinant of these equations are solved [29] to get 
dispersion relations. Since the determinant of the coefficients of the Fourier analyzed equations cannot be solved 
analytically, wc solve it with the help of software Mathematica for special cases and wave numbers are obtained 
in terms of angular frequency and lapse function. We use wave number to calculate phase and group velocities, 
refractive index and change in refractive index w.r.t. angular frequency. The graphs of each of these quantities 
are given and discussed. 

The paper has been organized as follows. In the next section, we shall review the general results for the 
3+1 formalism. The GRMHD criuations in the Schwarzschild black hole magnetosphere for the cold plasma 
are also given. In section 3, we restrict the problem only to the non-rotating background. Section 4 is devoted 
to the rotating non-magnetized plasma and section 5 is furnished with the dispersion relations for the rotating 
magnetized plasma. In the last section, we shall summarize and discuss the results. 

2 3+1 Spacetime Modelling 

In the 3 + 1 formalism, the line element of the spacetime can be written as [19] 

ds^ = -a^dt^ + 7y {dx* + 13' dt) {dx^ + P^dt) , (2.1) 
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where a (lapse function), /3' (shift vector) and (spatial metric) are functions of the coordinates t, . A nat- 
ural observer, associated with this spacetime called fiducial observer (FIDO), has four velocity n perpendicular 
to the hypersurfaces of constant time t and is given by 



Throughout the paper, geometrized units with G = c = 1 will be used. Vectors and tensors living in four- 
dimensional spacetime will be denoted by boldface italic letters such as FIDO's four-velocity n. The three- 
dimensional tensors are distinguished from vectors by a dyad over the letter, such as three-dimensional metric 
All vector analysis notations such as gradient, curl and vector product will be those of the three dimensional 
absolute space with three metric The determinant of three metric is denoted as g: g = rfet|7ij|. Latin letters 
i,j, fc, ... represent indices in absolute space and run from 1 to 3. 
The Schwarzschild line element in Rindler coordinates [9] 

ds^ = -a'^{z)df + dx^ + dy^ + dz^ (2. 3) 

is the best approximation of this spacetime in Cartesian coordinates. Since the Schwarzschild black hole is 
non-rotating, hence the shift vector /? = 0. This model is an analogue of the Schwarzschild spacetime with z, 
X and y as radial r, axial <f) and poloidal 6 directions respectively, a is the lapse function which depends on r 
in the Shwarzschild metric, whereas it depends on z in this analogue. In Rindler coordinates, a = -^^i where 
Th is the radius of the Schwarzschild black hole. This function vanishes at the horizon which we can place at 
z = Q and it increases monotonically as z increases from to oo. The associated FIDO's four-velocity vector 
reduces to n = The constant time surfaces have spatial geometry ds^ = dx^ + dy^ + dz"^ which gives a 

good approximation to the Schwarzschild geometry near the horizon. 



2.1 The 3+1 Split: Laws of Physics 

The real power of 3-1-1 viewpoint lies in its strong similarity to a physicists with non-general relativistic ex- 
perience. The 3+1 laws of electrodynamics arc formulated with an eye towards emphasizing this similarity. 
However, there are real unavoidable and crucial new features of these familiar laws that enter electrodynamics 
in the vicinity of a black hole. The most important of these is the issue of two times, the universal time t and 
the FIDO's time t. Both times arc related to each other by a = For the Schwarzschild black hole, the 
Maxwell equations in 3-1-1 formalism take the form [9] 

V.B = 0, (2. 4) 

— = -Vx(aE), (2.5) 
V.E = 47rpe, (2. 6) 



'dt 



V X (aB) - 47rQ:j, (2. 7) 



where B and E are the FIDO measured magnetic and electric fields, j and pe are electric current and electric 
charge density respectively. Under the perfect MHD condition, i.e., 

E-l- V X B = 0, (2. 8) 

where V is the velocity of the fluid measured by the FIDO, there can be no electric field in fluid's rest frame. 
For perfect MHD, the equations of evolution of magnetic field [18] become 

^ + -(B.V)(/3 - aV) + [e + -V.(aV)l B = 0, (2. 9) 

Ut a \ a } 

where 
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The expansion rate of FIDO's four- velocity is 



Here dot indicates partial derivative with respect to time t. The FIDO's measured rate of change of any three 
dimensional vector in absolute space (orthogonal to n) is 

dr dt ^ ^ 

Local conservation law of rest-mass according to FIDO is given by [18] 

''^°+P07^V.^ + ^|f + V.(aV-/3)|=0, 



Dt Dt a [2g 

or 



n ■ Pcn'V.— + po <^ ^ + -V.(aV) ^ = (2. 10) 

Dt Dt {a J 

with po as the rest mass density. FIDO measured law of force balance equation can be written as 

= {ai - l^f3j\iV^ - i£aij)V^^ -p\i + i-(V X B)iV.(V x B) - ^(aB)^ 

+ -^(aBi)|j-B'' - -^(B X {V X [V X (aV X B) - (B.V)/3] + (V x B).V/3})i (2. 11) 

with J), the pressure. Notice that the acceleration, a = is a function of time lapse in 3+1 formalism. 

For the Schwarzschild black hole with perfect MHD assumption given by Eq.(2.8), the equations of evolution 
of magnetic field (2.5) and (2.9), the mass conservation law (2.10) and the momentum conservation law (2.11) 
take the following forms respectively 

r)B 

_ = -V X (aV X B), (2. 12) 

— + (aV.V)B - (B.V)(aV) + (V.(aV))B = 0, (2. 13) 

at 

+ {(aV).V}(poM) + P0M7'V.^ + poM7'V.(aV.V)V + po/x{V.(aV)} = 0, (2. 14) 

+P07'V, + (V.V)m| = -poM7'ai - + ^(V X B),V.(V x B) - ^(c^B)^ 

+ ^(ai?0,iB^' - ^[B X {V X (V X (aV X B))}],. (2. 15) 

Eqs.(2.12)-(2.15) are the perfect GRMHD equations for the Schwarzschild black hole. In the rest of the paper, 
we shall analyze these equations using perturbation and Fourier analysis procedures. 



2.2 Cold Plasma and GRMHD Equations 

Theoretical modelling of the moving plasmas neglects the thermal and gravity effects (such as the pressure forces 
acting on plasma may not be as important as electromagnetic and centrifugal forces). This plasma is called cold 
relativistic plasma. This is the simplest closed system and contains only the equations of conservation of mass 
and momentum. The highest moment of this system, i.e., the kinematic pressure dyad is taken to be zero. This 
model can be used in the study of small amplitude electromagnetic waves propagating in plasmas with phase 
velocity much larger than the thermal velocity of the particles. 
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In hydrodynaniic trcatiiiciit. the plasma is represented as a perfect fluid which has no viscosity and no heat 
conduction. In Rindler analogue of the Schwarzschild spacetime, we assume that the system of equations for 
perfect MHD is enclosed by the cold plasma. This state can be expressed by the following equation [19] 

IJL = — = constant. (2. 16) 

We note that the cold plasma has vanishing thermal pressure and vanishing thermal energy. Here p is the mass 
density of the fluid. 

Using Eq.(2.16), the GRMHD Eqs.(2.12)-(2.15) become 

— = V X (qV X B), (2. 17) 
FiTK 

— + (aV.V)B - (B.V)(aV) + BV.(aV) = 0, (2. 18) 

1^ + (aV.V)p + P7^V.^ + p72V.(aV.V)V + pV.(aV) = 0, (2. 19) 

{ + 1) + - h^^^] {hit + ^•^) - i^'^^ - h^^^) ^^^^'^ + ^'^^ 

liW X B).V.(V X B) - ^ + Mik^ - ^[B X {V X (V X (aV X B))}]. (2. 20) 

The perturbed flow in the magnetosphere shall be characterized by its velocity V, magnetic field B (as measured 
by the FIDO) and fluid density p. The first order perturbations in the above mentioned quantities are denoted 
by 5V, (5B and 5p. Consequently, the perturbed variables will take the following form: 

B = B° + Bb, V = V°+v, p = p° + pp, (2.21) 

where B°, V° and p° are unperturbed quantities. Since waves can propagate in ^-direction due to gravitation 
with respect to time t, the perturbed quantities depend on z and t. 



3 Non-Rotating Background 

In this section, we shall solve GRAfflD equations by taking non-rotating background. Notice that this makes 
no difference whether we take magnetic field zero or non-zero. The relative assumptions for this background 
are given. The perturbations and Fourier analyze techniques are used to reduce GRMHD equations to ordinary 
differential equations. The numerical solutions of dispersion relations are given in the form of graphs. 



3.1 Relative Assumptions 

For non-rotating background, the magnetosphere has the perturbed flow only along z-axis. The fluid four- 
velocity measured by FIDO is described by a spatial vector field lying along the ^-axis and is given by 

V = u{z)e^. (3. 1) 

The Lorentz factor 7 = , ^ which takes the form 

FIDO measured magnetic field is given by 

B = B(z)e^. (3. 3) 

The following notations for the perturbed quantities will be used 

, , , 
b = =h^{t,z)e^, 

V = = Vz(t, -2)ez, 

p ^ ^-^=p{t,z). (3.4) 
P 
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We also assume that the perturbations have harmonic space and time dependence, i.e., 

/5, b, V - e-*('^*-'=^\ (3. 5) 

Thus the perturbed variables can be expressed as follows: 

b,(i,^)=C3e-'('^*-'=^), (3.6) 

where Ci, C2, C3 are arbitrary constants. 
3.2 Perturbation Equations 

We use linear perturbation to write down the GRMHD Eqs.(2.17)-(2.20) by using Eq.(2.21) 

= V X (av X B) + V X (aV X ^B), (3. 7) 

V.(^B) = 0, (3. 8) 

V.V)5p + P7'V.(-^ + V.V)v - ^(V.V)p + p(V.v) 

a ot am p 

= -2p72(V.v)(V.V);n7- p72(V.VV).v + p(v.VZTOi), (3. 9) 



= -72{<5p + 2p72(V.v)}a, + —{{aB,),^ - {aBj)^i}dB^ 



1 

Ana 

-fry'^iviV^ + v'Vi)VkjV'' - -f^{6pV^ + 2p-f^{\'.v)VJ + pv^jV^j 

--f%{6pV^ + WCV.v)^^' + pv^}Vj,kV''. (3. 10) 

It is to be noted that the law of conservation of rest-mass [18] in three dimensional hypersurface 

apju = constant 

is used to obtain Eq.(3.9). This will also be used in deriving the component form of these equations. 
The component form of Eqs.(3.7)-(3.10) is given as follows: 

-§ = 0, (3.11) 
a at 

bz,, = 0, (3. 12) 

-f + up^z + I'u-^ + (1 + ^'u')v.,. = (1 - 27V)(1 + j'u')-Vz, (3. 13) 

a ot a ot u 

7^(1 + 72^2)!^ + ^2^^ ^ ^^u^)uv,,, = -n^{a, + u{l + ^^u^)u'} 
--y^Wil + 7^M^)(1 + 47V) + 2u'y^a,)}v,. (3. 14) 

When the above equations are Fourier analyzed by substituting values from Eq.(3.6), these take the form 

-— C3 = 0, (3.15) 

a 

Lkcs = 0, (3. 16) 



ci + tku] +C2< + )i'k-(l-2Yn )il+Tu^) l u} =0, (3. 17) 

\ a J I u a ) 

ci7^{a^ + uu'{l + 7^w^)} + C2[7^(l + 7^"^) + ''Uk^ + 

+7^{u'(l + 7\^)(1 + 47^^) + 2u7^a4] = 0. (3. 18) 
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Eqs.(3.15) and (3.16) show that C3 is zero which imphes that there are no perturbations occurring in magnetic 
field of the fluid. Thus wc are left with Eqs.(3.17) and (3.18). It is mentioned here that for the non-magnetized 
plasma, we also obtain the same two equations. 

3.3 Numerical Solutions 

We consider time lapse a = z and p as a constant. Using mass conservation law in three dimensions, we obtain 
the value of u = . When we use these values, the determinant of the coefficients of constants ci and co 

in Eqs.(3.17) and (3.18) give us a complex number. Equating it to zero, wc obtain three vahics for fc, two from 
the real part and one from the imaginary part. Using the value of k, one can calculate the following quantities: 

1. Phase Velocity {vp): The velocity with which the carrier wave in a modulated signal moves and can be 
calculated by ^. 

2. Refractive Index (n): A property of a material that changes the speed of light, computed as the ratio 
of the speed of light in a vacuum to the speed of light through the material. Refractive index is the inverse 

of the phase velocity. 

3. Change in refractive index with respect to angulcir frequency (3^): This term helps to find out 

whether the dispersion is normal or not. 

4. Group Velocity {vg): This is the speed of transmission of information and/or energy in a wave packet 
and can be calculated by the formula — . 

Wc would like to mention here that only longitudinal oscillations are calculated in this case as we are using 
the longitudinal part of the momentum equation. Thus we would discuss the longitudinal waves propagating 
parallel to the magnetic field B in the remaining part of this section. 

The dispersion equation from real part of the matrix is of the form A(z)k^ + B{z)kLU + C{z)uj^ + D(z) — 0. 
This leads to two values of k in terms of z and angular frequency shown in the Figures 1 and 2. These are 
conjugate roots of quadratic equation. The real wave numbers are shown in the graph lying in the region 
3.5 X lO"'^ <z< 1.41. 

We see from the Figure 1 that angular frequency lo and wave number k are directly proportional to each other, 
i.e., the change in angular velocity gives the respective change in wave number. Also, the wave number decreases 
as z increases. This means that waves lose energy as wc go away from the black hole horizon. At the horizon, the 
wave number becomes very large and the waves vanish due to the effect of gravity. There is an abrupt increment 
in phase velocity when the angular frequency is near to zero. The region 3.5 x 10~^ < z < 0.5, 1.5 < a; < 10 has 
refractive index greater than one and its variation with respect to omega greater than zero. This corresponds 
to the region of normal dispersion [30]. The group and wave velocities go parallel and are increasing as we go 
away from the horizon. 

The Figure 2 shows that the waves are growing with the increase in angular velocity, but damping occurs 
with as z increases. For < z < 3.5 X 10"^ the wave number goes to infinity and hence there is no wave at the 
horizon. This shows that the event horizon is the interface at which they are formed and as we go away from 
the event horizon, the intensity of these waves decay. Since the refractive index is greater than one and ^ is 
greater than zero, the whole region shows normal dispersion. 

The dispersion equation obtained from the imaginary part of the matrix is of the form A{z)k + B{z)u) = 0. 
Hence we are with only one value of k given by the Figure 3. 

The Figure 3 shows that the waves gain energy with the increase in angular velocity, but lose as z increases. 
The wave number goes to infinity and hence there is no wave for < 2: < 3.5 x 10~^. The group and phase 
velocities are showing the same behavior i.e. they are increasing with an increase in z. The refractive index is 
not greater than one and hence the dispersion is not normal. 
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Surface Plot of Wave Number 




Figure 1: The waves decrease as they move away from the event horizon. The phase and the group velocities 
are equal except for the waves with less angular frequencies. The region shows normal dispersion. 
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Surface Plot of Wave Number 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 2: The waves damp as they go away from the event horizon. The phase and the group velocities are 
equal except for the waves with negligible angular frequencies. The region shows normal dispersion. 
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Surface Plot of Wave Number 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 





Figure 3: The waves damp as they move away from the event horizon. The phase and the group velocities are 
same. The dispersion is not normaL 
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4 Rotating Non-Magnetized Background 

When we consider non-magnetized cold plasma in rotating background, i.e., B = 0, the GRMHD Eqs.(2.17) 
and (2.18) vanish. Eqs.(2.19) and (2.20) change into general relativistic hydro dynamical equations given below: 

^ + (qV.V)p + P7^V.^ + frr^Y.{a\.V)y + pV.(aV) = 0, (4. 1) 



{P7%j + P^^Vj) ( -- + V.V ) = -pjW (4. 2) 



ld_ 

a dt 

For rotating background, the fluid four-velocity measured by FIDO can be described in x^r-plane as follows 

V = V{z)e^ + u{z)e^ (4. 3) 

and the Lorentz factor 7 takes the form 

Here we shall use the following notations for the perturbed quantities 

p = —= p{t, z) 

V = 5Y = Vx(i, z)e^ + v^{t, z)e^ (4. 5) 

and the perturbations have harmonic space and time dependence, i.e., 

p, V ~ e-'('^*-'=^\ (4. 6) 

The perturbed variables take the form 

p(i,^) = Cie-'("'-'=^), 
Vz(i,2) =C2e-'('"*-'=^), 

Vx(i,2)=C3e-'("*-'=^). (4.7) 

4.1 Perturbation Equations 

Introducing linear perturbations in Eqs.(4.1) and (4.2), it follows that 

i| + V.v) 6p + p,^V. (1| + V.v) V - ^(V.V)p + p(V.v) 
= -2p72(V.v)(V.V);n7 - p^'^{Y.VY).w + p{v.Vlnu), (4. 8) 

irr'Sij + p^%Vj)^^ + P7\jV^ + P7%vj,kVW'' = -7'{5p + 2p'y\V.v)}ai 
-p-f\viV^ + v^Vi)Vk,jV'' - ^^{6pV^ + 2pj^{\.v)V^ + pv^jVij 

--f%{SpV^ + 4p72(V.v)y^' + pv^jVj^kV''. (4. 9) 



11 



The component form of Eqs.(4.8) and (4.9) gives the following three equations 

1 dp Itr f ^ '^^x \ 2 1 (^''^z 2 2\ 

+Up,z + J V ^+UVx,z +7 ^ + (1 +7 U )Vz,z 

a at \a ot J a at 

= -7^u{(l + 2-f^V^)V' + 2-f^uVu}v^ + |(1 - 2^'^u^){l + 7^u^)^ - \ v^, 

2/1 , 2t^2n f ^ dVa: \ 4 1 St;^ 

= -P7^u{(l + j'^V^)V' + ^^uVu'} - + 4^^V^)uu' + AVV'{1 + ^^V^)}v^ 

+ 27\2)(i + 27^2) - j'^V^jV + 27^(1 + 2j\'^)uVu']v„ 

7^(1 +7V) (i^+uv^^ +^\V (^^^+uv^„ 
\a at J \a at 

= -P7^{a^ + (1 + -i'^u'^)uu' + -i^u^YV'} 
-[7*{uV(l + 472y2) + 2V{a, + (1 + 2-i'u^)uv!)}\v^ 
-^\v!{\ + 7^«^)(1 + 47^w^) + 2w7^{a, + (1 + 2^'^u^)VV'}\v^ 
Using Eq.(4.7), the Fourier analyzed form of Eqs.(4.10)-(4.12) become 



Cl 



■ C2 



-— 7'm + (1 + ^''u^yk - (1 - 272i,2)(i + ^2^2^^ ^ 27VyF' 



+C37 



to; 



y + LkuV + 72m{(1 + 27^^)1/' + ^-i^uVv!} 



+{(1 + 272u2)(i + 272^2) _ ^2y2|-^/ ^ 272(1 + 2^'^u^)uVv!\ 

+C3[(-^ + t/fcw) 72(1 + 72^2-) ^ ^4^|(j ^ 472f2)mm' + 4Fy'(l + 7^^)}] = 0, 

cn^Oz + (1 + l^u^W + 72M2y]/'} + C2[7^(l + l^u^) (-^ + '-^^^) 
+72{u'(l + 72-Lt2)(l + 472^2) + 2u72(a^ + (1 + 2'fu^)VV')\\ 



+C37 



(-^ + mfc) ul/ + u^V(\ + 472^2) + 2y{a^ + W(l + 272^2)} 



0. 



(4. 10) 



(4. 11) 



(4. 12) 



(4. 13) 



(4. 14) 



(4. 15) 



The determinant of the coefficients of Ci, C2 and C3 gives a complex number. Equating it to zero yields a 
complex dispersion relation in k. 



4.2 Numerical Solutions 

We consider time lapse a — z and assume that V = u for the sake of convenience. Using mass conservation 
law in three dimensions (with mass density as a constant quantity) we obtain u = ■ ^""^ ^'^^^ section we are 

having modes when B = 0. The complex dispersion relation obtained from the determinant of coefficients leads 
to two dispersion equations. The real part gives an equation of the type A{z)}p' + B{z)ku) + C{z)io^ + D{z) = 
leading to two dispersion relations shown in the Figures 4 and 5. The imaginary part shows a cubic equation 
A{z)k^ + B{z)k^LV + C{z)koj'^ + D{z)oj'^ = 0. This equation has one real and two complex conjugate solutions. 
The Figure 6 shows the real value of k obtained from this equation. 

The Figure 4 indicates that the region for which real waves exist lies in 1.5 x 10"^ < z< 10, 0.045 < w < 10. 
It is clear from the Figure that the region < w < 0.045 contains evanescent waves. We note that at the horizon 
and nearby values of z i.e. in the region < < 1.5 x 10~^, the wave number is infinite and hence no wave 
exists there. The wave number is decreasing as we go away from the horizon and hence damping occurs. Also 



n > 1 



dn 
' dui 



> and hence the region is of normal dispersion. The phase velocity is less than group velocity of 
waves except at some points where the refractive index n is positive but is less than one, as shown in the graph. 
This means that the dispersion is not normal but anomalous at these points [31], normal otherwise. At a; = 
and very small angular frequencies, the phase velocity is high but the group velocity become infinite which is 
not physical. 
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Surface Plot 



of Wave Number 




Figure 4: The waves damp when they move away from the event horizon. The phase velocity is less than the 
group velocity in most of the region. Most of the region shows normal dispersion. 
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Surface Plot of Wave Number 




The Figure 5 contains the infinite wave number in the region < z < 1.5 x 10^^, that is no wave present 
there. The wave number possesses real values in the region 1.5 x 10~® < z < 10, 0.05 < lu < 10. On increasing 
the value of z, the waves arc; damping but are growing with the increase of the angular frequency. In this 
region, refractive index n < 1 and < which implies that the region is not of normal dispersion. At 
1.5 X 10~^ < ^; < 10, < a; < 4, the wave number becomes imaginary at several points leading to the result 
that evanescent waves present there. 

In the Figure 6, real waves exist in the region 1.5 x 10^^ < z < 10. The refractive index is greater than 
one, ^ is negative and thus the region is not of normal dispersion. No wave is present on the event horizon 
because the wave number is infinite there. As we go away from the horizon, the wave number decreases. On 
increasing the angular frequency, we find growing waves. Also, group velocity is greater than phase velocity 
which indicates that the region is of anomalous dispersion. 

5 Rotating Magnetized Background 

The GRMHD equations for magnetized cold plasma in rotating background will remain the same as given by 
Eqs.(2.17)-(2.20). Here the four-velocity of the fluid measured by FIDO can be described in the same way as 
given in the previous section. The rotating magnetic field can be expressed in the xz-plane as 

B = S[A(z)ex+e,], (5.1) 

where A, u and V are related to each other by 

V=— + Xu, (5. 2) 

a 

Vp as an integration constant. The following notations for the perturbed magnetic field will be used in addition 
to the notations given by Eq.(4.5). 

b = ^ = bx(t,z)ex + bz(t, z)ez. (5.3) 

Since the perturbations have harmonic space and time dependence, i.e., b ^ Q-i-i'^t-kz) ^ ^^i^ perturbed variable 
can be expressed as follows 

bx(t,^) = C4e-'('^*-*=^), 

b,(t,z) = C5e-'('^*-'=^). (5.4) 

5.1 Perturbation Equations 

Using Eq.(2.21) in Eqs.(2.17)-(2.20), we have 

^^51 = V X (av X B) + V X (aV X ^B), (5. 5) 

V.(^B) = 0, (5. 6) 



-757 + ^-^ H-" + P^^- -7^7 + h - — (V-V)p + /9(V.v) 

a at J \a ot J p 

= ~2pj^(V.v){V.V)lnj - p72(V.VV).v + p{v.Wlnu), (5. 7) 

+Pl^v^,JV^ + PJ^V.Vj^kVW" - -^{{adBij^j - (aSB/),,}B^ = -^^Sp + 2pj^{V.v)}ai 

47ra 

+ J-{(^aB,),j - {aBj)^i\5B^ - p^\v,V^ + t^^VO^fc.j^' " l^i^V + Ipi'iy .-f)V' + tyv'}Vr., 
47ra 

--i^Vi{bpV^ + 4p72(V.v)y^' + pv^}Vj^^V^. (5. 
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Surface Plot of Wave Number 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 6: The waves damp as they move away from the event horizon. The group velocity is greater than the 
phase velocity. The region is of anomalous dispersion. 
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The component form of the Eqs.(5.5)-(5.8) can be written, after a tedious algebra, as follows 

1 db 
a at 

-^+uhz,z = 0, 
a at 



Kz = 0, 

-— + +7 ^--^ + 7 ^~~g^ + (1 + 7 u )Vz,z +7 uVvoc,z 



1 

a dt 



a dt 



(5. 9) 

(5. 10) 
(5. 11) 



+ 27^2)^' + 272uyu>^ + [(1 - 2^^u'){l + ^"^u^)- - 2^^u''VV'\v„ 

u 



1 1 dVr 



All \ a dt 



fTf uV ■ 



AB2\ Idv 



+ pj^V + 



An 



UV 



47r J a dt 
52 



+ |p7'(l+7V2) 
bx{oe'{l — u^) — auu'} 



B^ 
An 



+PPYu{{1 + YV^)V' + YuVu'} + 



P7'^u{(i + AYv^)uu' + 4(1 + Yv )yy'} + 



B'^ua'' 
Ana 



(5. 12) 



+ 



B^u. 



P72[{(1 + 272u2)(l + 27^2) - 72y2}y' + 27^(1 + 27V)uyu'] + ^(Aa)' 



P7^(l +7V) + 



A^B^ 1 1 dv 



, A ir XB^ 

+ ( P7 + -^J'fJ'Vx 



An I a dt 

XB^ 



pj uV — 



XB^\ 1 dvx 



An J a dt 
B^ 



Ana 

^7^(1 + 7V) - 



Vz=0, 



(5. 13) 



4-77 



■ UV, 



^ (1 - u'^)bx,z + -^{{aXy + a'X - uX{ua' + u'a)}bx 
An Ana 



^pY\az + w{(l + 7 w^)"' + ruVV'}\ 



+ 
+ 



P7''{w V(l + A-fV^) + 2V{az + uu{\ + 27^^))} 



47ra 



P7^{u'(l + 7^u^)(l + 47V) + 2u7^{(l + 27V)yF' + a J} 



Ana 



iaxy 



Vz = 0. 



(5. 14) 
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When the above equations are Fourier analyzed, they take the following form 

C3(q;' + ika) — C2{{aXy + LkaX} + C5{aVy — C4{(aw)' — llv + tkau} = 0, (5. 15) 

C5 (-— + ikuj =0, (5. 16) 

C5ik = 0, (5. 17) 
ci (-if^ + ,ku) + C2 (--7'« + ik{l + 7'«2) - (1 - 272«2)(i + ^2^2)^ ^ 27Vyy' I 

V Of / I W J 

+C37^ (-^ +tfcu) y + u{(l + 27V^)y' + 272uyu'} =0, (5.18) 

01/0724(1+721/2)^' + ^^uVu'} - ^C4{(1 - u^)Lk + —(1 - u2) - W} 

47r a 



+C2 



- ( P'y^uV - ^) — + iku{ fTf^uV + ^ ) + P7'{(1 + 27'w')(l + 272^2) _ ^V2}y' 



47r / Q 



47r 



+2^7^(1 + 272u2)ul/u' + -^(Aa)' 
^ ^ 47ra ^ ^ 



+ C3 



|p7'(l + l^V^) -^\+ P7^w{(l + 47^2) W + 4(1 + 72y2)yy'} _ 

y Ait \ Ana 

CiP7'[a. + u{{l + -f^u^)u' + j^VuV'}] 



0, (5. 19) 



+C2 



-<^P7ll+7V) + 



An I a 



\2 r2 



^7^(1 +72^2) - 



47r 



XB^u, 



+ { n^u'il + 72w2)(i + 4^2^2) ^ 2u72((i + 272^2)^^' + a,)} - ^^(aA)' 



+C3 



— ( p7^uy — — ^ ) — + ( P7^uy + 1 tfcu 



47r / a 



47r 



+ <^ p7'*Ky'(l + AYV^) + 2V{a, + uu'{l + 2Yu^))} + 



XB^a'u 
Ana 



£2 
+— C4 

4-77 



A(l - u^)ik + A-(l - ^2) _ Xuu' + ^ 
a a 



= 0. 



(5. 20) 



Equation (5.17) implies that C5 is zero which gives that bz is zero. 



5.2 Numerical Solutions 

We consider the same circumstances as wc; have considered in the rotating non-magnetized plasma and assume 



that u 



V. Substituting this value in Eq.(5.2) with Vp = 1, it follows that A = 1 



vW2 



which 



shows that the magnetic field diverges near the horizon (as discussed in [25]). For the sake of convenience, we 
take B = 



so that ^ = 2. Making use of C5 = (Eq.(5.17)), we get a 4 x 4 matrix of the coefficients 
of constants. The complex dispersion relation has two parts. The real part gives a dispersion relation of the 
type ^(z)fc4 + B{z)k^uj + C{z)uj^ + D{z)uj'^ + E(z)k:^ + F{z)k^uj'^ + G{z)kuj + H(z)kuj^ + L{z) = whereas the 
imaginary part gives a relation of the form A{z)k^ + B{z)k'^LU + C{z)uj + D{z)uj-^ + E{z)kuj'^ + F{z)k = 0. The 
first equation gives four dispersion relations out of which two are real and interesting while the other two roots 
are complex conjugates of each other. The second equation gives three dispersion relations. The respective 
wave numbers are shown in the Figures 9-11. These relations give different types of modes when B > and 
wave number is in arbitrary direction to B. The two dispersion relations obtained from the real part are shown 
in the Figures 7 and 8. 

The Figure 7 shows that at the horizon the wave number is infinite and hence no wave is present there, k 
has real values in the region 1 < 2; < 10, 10^'' < w < 10. No wave is present in the region < 2; < 1, evanescent 
waves are present in < a> < 10~^. All the quantities, the wave number, wave velocity, group velocity and 
refractive index are negative. This leads to the existence of metamaterial in the region [32]. The wave number 
is decreasing with an increase in angular frequency and increasing with an increase in time lapse. 
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Surface Plot of Wave Number 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 7: The waves grow as they move away from the event horizon. The medium admits the properties of 
metamaterials. 
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Surface Plot of Wave Number 




The Figure 8 shows that the wave number increases with the increase in angular velocity and decreases 
with the increase of lapse function. The horizon possesses no wave due to infinite wave number. The region 
< 2; < 1 has no wave and < uj < 5 x 10^** is region of evanescent waves. The value of refractive index is less 
than one and hence it is not a normal dispersion. 

The dispersion relation obtained from the imaginary part of the matrix determinant are shown in the Figures 
9, 10 and 11. 

The Figure 9 shows that the wave number is infinite in < z < 0.5. The wave number is positive in the 
region 0.5 < 2; < 1.6, it decreases abruptly and then it smoothly rises but is negative. The phase and group 
velocities are negative. The refractive index is negative in the domain 1.6 < 2; < 10 so the region possesses the 
properties of metamaterials. 

In the Figure 10, real waves exist in the region 0.5 < 2 < 10. No wave is present on the event horizon as 
well as in the region < 2; < 0.5 because the wave number goes to infinity there. The waves are increasing in 
0.5 < 2 < 0.7 then there is a sudden decrease. The region 0.5 < 2 < 0.7, < w < 4 contains evanescent waves. 
The refractive index is less than one in the whole domain, hence the waves are not normally disperse. 

The Figure 11 shows that real waves exist in the region 0.5<2<10, 0.5<w<10. The phase and group 
velocities are same in this domain. The refractive index is less than one and the dispersion is not normal. No 
wave is present on the event horizon and the region near the horizon because the wave number is infinite there. 
As we go away from the horizon, the wave number decreases. On increasing the angular frequency, we find 
growing waves. 

6 Conclusion 

We have reformulated the GRMHD equations for the Schwarzschild black hole magnetosphere to take account of 
gravitational effects due to event horizon. We discuss one dimensional perturbations in perfect MHD condition. 
This has been explored for the cold plasma case. These equations are written in component form to make them 
easier for the Fourier analysis. The equations are Fourier analyzed by using the assumption of plane wave. The 
determinant of the coefficients is solved numerically for the following backgrounds in the magnetosphere. 

1. Non- Rotating Background (either non-magnetized or magnetized) 

2. Rotating Non-Magnetized Background 

3. Rotating Magnetized Background. 

The non-rotating background shows pure Schwarzschild geometry outside the event horizon whereas the rotat- 
ing background demonstrates the restricted Kerr geometry in the vicinity of the event horizon which admits a 
variable lapse function with negligible rotation. We have taken Schwarzschild black hole with rotating back- 
ground due to the well-known fact that the Schwarzschild black hole is the minimum configuration of the Kerr 
black hole. 

In the case of non-rotating background, we have found that the magnetospheric fluid is normally disperse. 
The wave number becomes indefinite at the horizon and for the values of z near the horizon and hence no wave 
is present there. This shows that no signal can pass the event horizon or near to it. The third subcase shows a 
mode which is not normally dispersive. 

In the case of rotating non-magnetized background we have found one case (Figure 4) in which the dispersion 
is normal. The regions with zero angular frequency or very less angular frequency are evanescent in all the three 
cases. The wave number becomes infinite at event horizon and in the neighborhood of the event horizon and 
no wave exists there. 

In rotating magnetized background, we obtain region with negative wave number, phase velocity, group 
velocity and refractive index (Figures 7 and 9). This fact is deduced by the analysis of two of the dispersion 
relations obtained in section 5. The refractive index is negative in both the cases which implies that the region 
possesses all the qualities of left-handed metamaterials. The other three dispersion relations give that the region 
under discussion is of non-normal dispersion which restrict the real signals pass through that region. 

It is worth mentioning that the wave number becomes infinite at the event horizon and no wave exists 
there. This supports the well-known point of view that no information can be extracted from a black hole. It 
is interesting to mention that this work extends the results given by Mackay et al. in [33] according to which 
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Surface Plot 



of Wave Number 




Figure 9: The waves damp and then grow as they move away from the event horizon. The medium admits the 
properties of metamaterials. 
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Surface Plot of Wave Number 




Figure 10: The waves suddenly grow and then damp on going away from the event horizon. The dispersion 
not normal. 
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Surface Plot of Wave Number 





Figure 11: The waves damp when they move away from the event horizon. The phase and the group velocities 
are equal. The dispersion is not normal. 
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rotation of a black hole is required for negative phase velocity propagation. It is also observed that the waves 
of less angular velocity are evanescent. Our numerical results indicate that negative phase velocity propagates 
in the rotating background whether the black hole is rotating or non-rotating. 

We know that the MHD waves in cold plasma are non-dispersive. However, the dispersion is noted in 
the above figures. This factor comes due to the formalism used and the equations which provides different 
equations from the usual MHD equations. Since the 3-1-1 split of GR is used in a preliminary investigation of 
waves propagating in a plasma influenced by the gravitational field. Thus internal gravity waves which interrupt 
the MHD waves imply the cases of dispersion in each of the hypersurface. We would like to mention that graphs 
of the waves are given in a particular hypersurface where time is constant, not in the whole Schwarzschild 
background and consequently this is justified locally, not globally. Finally, it is mentioned that some of the 
Figures have patches missing which is due to the existence of complex numbers there. Mathematica cannot plot 
complex numbers with real numbers. These complex numbers are shown by gaps in the Figures. 

This analysis has been done for the cold plasma state. It would be interesting to extend this analysis for 
the isothermal state of plasma. This would provide the effects of pressure in the current work. Currently, it is 
in progress. 
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